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Abstract 

We study the current and charge distribution of a two dimensional electron gas under strong 
perpendicular magnetic fields within the linear response regime. We show within a self-consistent 
screening theory that incompressible strips with integer values of local Landau-level filling factor 
exist for finite intervals of the magnetic field strength B. Within an essentially local conductivity 
model, we find that the current density in these B intervals is confined to the incompressible 
strips of vanishing local longitudinal resistivity. This leads to vanishing longitudinal and exactly 
quantized Hall resistance, and to a nice agreement of the calculated Hall potential profiles with 
the measured ones. 
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I. INTRODUCTION 



Recent investigations [ll] with a scanning force microscope produced evidence, that the 
spatial distribution of the dissipative current, carried by a two-dimensional electron gas (2D 
EG) in a Hall bar under the conditions of the quantized Hall effect (QHE), is dominated 
by the existence of incompressible strips (ISs), which are expected to develop in the space 
regions where an integer number of Landau levels (LLs) is filled. |2] Subsequent theoretical 
work, ^ which calculated electron density and electrostatic potential within a self-consistent 
Thomas-Fermi-Poisson approximation, and current density and position-dependent electro- 
chemical potential from a local version of Ohm's law, confirmed the experimental evidence 
and found that the current flows in the ISs, in which the local longitudinal resistivity van- 
ishes. Due to the Thomas- Fermi approximation (TFA), incompressible strips were found for 
all values B of the magnetic field below a threshold Be, at which the center of the sample 
becomes incompressible, jsl] This leads to a vanishing total resistance for B < B^ to a poten- 
tial drop in several steps across ISs, which may coexist near a sample edge, and, for B < B^ 
to the absence of B intervals with linear variation of the Hall potential across the sample. 
These shortcomings are not in agreement with experiment 



of the model. 



l|] and require an improvement 



II. IMPROVEMENTS OF THE GUVEN-GERHARDTS MODEL 

We have recently improved on the approach of Giiven and Gerhardtsj^ in two essential 
respects^]. First, replacing the TFA by a Hartree-type approximation, which takes the finite 
extent of wavefunctions into account, we demonstrate that the existence of incompressible 

n 

strips is restricted to finite B intervals below -Be [4]. Second, we relax the assumption of 
a strictly local relation between conductivity and electron density, which leads to singular 
current density even along isolated lines of integer LL-filling factors, and we simulate non- 
local effects by a suitable averagingof the local conductivity tensor jj]. 

Following the previous work,j3|, |5| we calculate density profile and effective confinement 
potential of the 2D EG in thermal equilibrium self-consistently. Assuming that all charges, 
including a homogeneous positive background and induced charges on gates, defining the 
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Hall bar (in \x\ < d), are concentrated on the plane z = with translation invariance in y- 
direction, we calculate the potential V(x) from the electron density rieiix) solving Poisson's 

p n n 

equation under suitable boundary conditions. [2, |3|, |5|| Electrostatic self-consistency requires 
that the density of the spin-degenerate 2D EG is calculated from 

where / = ^yh/mujc is the magnetic length, Xq a center coordinate, /i* the electrochemical 
potential. The energy eigenvalues En{Xo) and eigenfunctions (f)n,Xo{x) should be calculated 
from the single-electron Schrodinger equation with the confinement potential V{x). In the 
previous work one assumed V{x) to vary slowly on the scale / and evaluated Eq.(0) in 
the TEA, En{Xn) fa E^, + V{Xo), |0„,Xo(a;)P ~ S{x - Xq), neglecting the extent of the 
wavefunctions. This yields flat ISs in the density profile with local filling factor 
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FIG. 1: (color online) (a) Calculated density profiles for the left half of a symmetric Hall bar. 
While the TFA yields incompressible strips with ^{x) = 2 and iy{x) = 4, only those with i/(x) = 4 
survive in the QHA. (b) Local filling factors (gray scale plot) versus position x and cyclotron energy 
ilc = ISs (light) exist in the TFA for all B with Q^c/E^ < 1, in the QHA, however, only for 
1.0 > njE^ > 0.65 and 0.5 > njE^ > 0.4. d = 3.1^m, average electron density 3 • lO^cm-^ 
temperature ksT/E^ = 0.02, E% = n%{Q)'Kh^ /m, where (0) = nei(O) at T = 0, = 0. 

v{x) = 2ttP nci{x) = 2 for all magnetic-field values B < i?2 (= Be), with z/(a;) = 4 for 
all B < -B4 (= B2/2), etc. This is, however, an artifact of the TEA. If one takes the 
width AvK of the wavefunctions into account, e.g. in the full Hartree approximation, one 
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finds that the ISs calculated in the TEA are smeared out if their width becomes smaller 
than Xw- This is already seen for the "quasi-Hartree approximation" (QHA), which uses 
again £'„(Xo) ~ + V{Xq) but replaces the (pn,Xo{^) by the Landau wavefunctions of the 
unconfined 2D EG, as is demonstrated in Eig^ and in Eig^, which shows within the QHA 
no ISs for = O.S^a <B < 0.6552. 

The stationary state with an applied dissipative current / was described Q by a local 
Ohm's law E(r) = /3(r)j(r), with p(r) = [(3"(?T,eir)]~-'^, which allows us to calculate current 
density j(r) and driving field E(r) = V/i*(r)/e with the constant components = and 
Ey{x) = Ey, and with jy{x) = Ey/pi{x), and Ex{x) = EyPni^x)/ pi{x), with pi{x) the 
longitudinal and Ph{x) the Hall resistivity. This yields the Hall resistance Rh = Vh/I and 
the longitudinal resistance Ri = 2dE^/I: 



The feedback of the applied current on the (measurable) change of the electrostatic 
potentialllj was considered by requiring the stationary state to be in local equilibrium. 
Here we consider only a linear response situation, in which the current-induced change of 
the electrostatic potential equals that of the electrochemical potential. To avoid spurious 
singularities of the integrals in Eq.(0) along isolated lines of vanishing pi{x), i.e. inte- 
ger which are an artifact of the strictly local model, we average the conductivity 
tensor, over a length A of the order of the mean electron distance (Eermi wavelength): 



Starting with the self-consistent Born approximation^] for the conductivity tensor of a 
homogeneous 2D EG, (3"(z/), replacing the filling factor u by the local z/(x) calculated within 
the TEA, and performing the spatial average, we find strips with vanishing denominators in 
Eq.(0) only for ISs which are wider than 2A and survive also in the QHA. If such wide ISs 
with z/(x) = 2k exist in a certain 5-interval, they carry all the current (cf. Eig|2K) and lead, 
for T — i> 0, to exactly quantized resistance values, Rh = h/{e'^2k), Ri = 0. Eor all other 
5- values the current is spread over the sample, and Eq.Q yields reasonable average values 
for Rh and Ri, as is demonstrated in Eig|2b- 





a{x) = 2^ 



J_,d^a{x + o. 
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FIG. 2: (color online) (a) Current distribution for selected B values. Only for il^c/Ep = 0.99 (solid 
line) the current is confined to the incompressible strips, which survived the averaging. A = 20nm. 
(b) Gray scale plot of the averaged i/(x) ^ (h/e^) (Th{x) versus B and together with Hall (solid) 
and longitudinal (broken light line) resistance; d = 7.5/im, A = 20nm, ksT/E^p = 0.02. 

III. SUMMARY 



In conclusion, the spatial averaging of the conductivity tensor, which simulates non- 
local effects and eliminates spurious singularities due to narrow ISs, which are an artifact 
of the TFA and do not survive in a Hartree-type approximation, represents an essential 
improvement of the theory and allows us to calculate reasonable results for the longitudinal 
{Ri) and the Hall {Rh) resistance as continuous functions of 5, with high-precision plateau 
values in the B intervals in which sufficiently wide ISs exist. Furthermore, we find three 
different types of the variation of the Hall potential across the sample (FigO)), in excellent 
agreement with the experiments. 3] 
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FIG. 3: (color online) (a) Plot of the scaled Hall potential ^{x)/<^>{d), <^{x) = 
versus B, together with RniB) (black) and Ri{B) (white line), d = 7.5/im, kT/E'^ = 0.0124. 
(b) $(x)/<I>((i) for selected values of B, vertically offset [cuts through (a)]. A stepwise variation is 
seen for B in a plateau of the QHE, a linear variation if B is far from a plateau, and a nonlinear 
variation in the center for B near a plateau edge. 
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